Simulating nanoscale dielectric response.

We introduce a constrained energy functional to describe dielectric response. We demonstrate that the local functional is a generalization of the long-ranged Marcus energy. Our reformulation is used to implement a cluster Monte Carlo algorithm for the simulation of dielectric media. The algorithm avoids solving the Poisson equation and remains efficient in the presence of spatial heterogeneity, nonlinearity, and scale dependent dielectric properties.